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Abstract: In a class of extensions of the minimal supersymmetric standard model with 
(B-L)/left-right symmetry that explains the neutrino masses, breaking R-parity symmetry is 
an essential and dynamical requirement for successful gauge symmetry breaking. Two con- 
sequences of these models are: (i) a new kind of R-parity breaking interaction that protects 
proton stability but adds new contributions to neutrinoless double beta decay and (ii) an 
upper bound on the extra gauge and parity symmetry breaking scale which is within the 
large hadron collider (LHC) energy range. We point out that an important prediction of 
such theories is a potentially large mixing between the right-handed charged lepton (e c ) and 
the superpartner of the right-handed gauge boson (W^), which leads to a brand new class 
of R-parity violating interactions of type ^ c V^e c and d c ^u c e c . We analyze the relevant con- 
straints on the sparticle mass spectrum and the LHC signatures for the case with smuon/stau 
NLSP and gravitino LSP. We note the "smoking gun" signals for such models to be lepton 
flavor/number violating processes: pp — > [i ± [i ± e + e~jj (or r ± r ± e + e~ jj) and pp — > ^e^bbjj 
(or r^e^bbjj) without significant missing energy. The predicted multi-lepton final states and 
the flavor structure make the model be distinguishable even in the early running of the LHC. 
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1 Introduction 

Supersymmetry (SUSY) is one of the popular and best motivated candidates for physics 
beyond the standard model (SM). It stabilizes the gauge hierarchy and provides a dark matter 
candidate in a natural manner. An intuitive requirement to stabilize the dark matter in 
MSSM is the existence of the R-parity symmetry under which all standard model particles 
are even and their superpartners are odd. The lightest supersymmetric partner field (LSP), 
e.g., either the neutralino or the gravitino, which is odd under R-parity is therefore suitable 
as the dark matter candidate. If R-parity is a global symmetry of the MSSM, it is logical 
to think of it as a remnant of some high scale physics. It will of course be interesting if the 
high scale physics is motivated by further reasons. A shortcoming of R-parity conserving 
MSSM is the zero neutrino mass. Understanding the origin neutrino masses then requires it 
to be part of a larger theory. An example of extension to the MSSM is to gauge the B — L 
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global symmetry, where anomaly freedom requires introducing a right-handed neutrino to 
each generation. The breaking of B — L symmetry gives Majorana masses to neutrinos. If the 
breaking is accomplished by Higgs fields with B — L = ±2, it not only helps to explain the 
small neutrino masses via the seesaw mechanism, but also leaves the R-parity as an unbroken 
symmetry at the level of the MSSM [1], thereby providing a stable dark matter candidate. 
Extending MSSM by a B — L symmetry therefore "kills two birds with one stone" . 

Two possible classes of models with B — L gauge symmetry are: (i) SU(2)l x U(1)y x 
U(1)b-l, and (ii) its left-right (LR) symmetric generalization based on SU{2)l x SU{2)r x 
U(1)b-l- Breaking B — L by two units in the second case is more appealing since it can 
explain the origin parity violation, and leads to a number of interesting phenomenological 
implications for LHC searches including the Wr boson as well for low energy weak processes. 
We discuss them in this paper. 

If the gauge symmetry is to be broken by a pair of Higgs superfields A c (l,3,+2) © 
A c (l, 3,-2) which are required to implement the seesaw mechanism and gauge anomaly can- 
cellation, two interesting results follow [2]. First, even though a priori the model is expected 
to have a remnant R-parity after symmetry breaking, in its minimal version, exactly the 
opposite happens, i.e., R-parity must be necessarily broken spontaneously in order for the 
full gauge symmetry to break down to the MSSM gauge group. If R-parity is exact, gauge 
symmetry cannot break [2]. If the model is extended to include singlets, there is a range of 
parameters where one can still have unbroken R-parity [3]. In the minimal model, however, 
R-parity breaking is mandatory. The right-handed (RH) sneutrino field, P c , has to pick up a 
vacuum expectation value (VEV) along with the neutral member of the B — L = 2 triplet, 
breaking the parity symmetry and contributing to the mass of the gauge bosons and gauginos 
associated with right-handed currents. Since v c is an R-odd particle, its VEV breaks R-parity. 
We call this class of models "dynamical R-parity breaking" models, since R-parity breaking 
is forced on the theory at the global minimum of the Hamiltonian. Other examples of models 
where R-parity breaking by v c vev are the minimal U(V)b-l extensions of the MSSM [4, 5]. 
In this note we will focus on the SUSYLR case. 

A consequence of dynamical R-parity breaking in minimal SUSYLR model is the pre- 
diction of an upper bound on the mass scale of the right-handed Wr boson, i.e., Jl% 8 < 
-^SUSy// 2 [6], which is in the range accessible at the LHC. Here Msusy is a generic soft 
SUSY breaking mass scale, / is the Yukawa coupling responsible for right-handed neutrino 
masses and has to be > 0.1. Similar relations are also found in SUSYLR models where v c 
vevs break left-right symmetry [7, 8]. 

Due to spontaneous R-parity breaking, neutrino masses arise not only from the usual 
type-I seesaw mechanism, but also via mixing with the neutralinos. Another consequence 
is that neutralino is no longer a stable particle and cannot therefore play the role of dark 
matter. However, if gravitino is the LSP, it can have an extremely long lifetime (> 10 26 sec) 
and play the role of dark matter [9]. Implications for such a dark matter particle have been 
studied extensively in connection with cosmic ray anomalies [10]. 

Since the scales of both superpartners and the new gauge interactions are predicted to lie 
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in the few TeV range, this theory could in principle be testable at the hadron colliders [11]. 
In this paper therefore, we study the genuine signals from dynamical R-parity breaking and 
discuss how it can be distinguished from usual R-parity breaking models at LHC. 

We point out that the most important consequence of dynamical R-parity breaking in 
SUSYLR models is a large mixing between the RH charged lepton and RH wino, i.e., the 
physical RH charged lepton after symmetry breaking is generically denoted as 

i c = e u £ c + e ew w+ + -- - , (i.i) 

where 0u,9ew ~ ^(1) for (u c ) ~ Mstjsy and the • • • represents the contributions of other 
Higgsino fields if the corresponding Higgses VEV's also violate parity. The physical charged 
lepton field contains a large RH wino component and in turn induces new R-parity violating 
terms of Kahler type. This is characteristic of dynamical R-parity breaking models [2, 7, 8] 
with the presence of gauged SU(2)r and it leads to new effects absent in usual R-parity 
violating MSSM or other models of spontaneous R-parity breaking, such as [12, 13]. In 
particular, it leads to effective R-parity violating interactions of the form /i c z/^e c , r c V^e c and 

d c \t c e c for all generations of quarks/squarks. We show that these kinds of vertices add new 
contributions to neutrinoless double beta decay and imply constraints on the parameters of 
the model. 

These new interactions bring rich phenomenology at the LHC. In the context of a realistic 
model based on left-right symmetry, we study the single production of a slepton NLSP from 
the RH neutrino decays, which is produced via an on-shell Wr boson resonance at LHC. The 
NLSP single production and decay yield multi-lepton final states of type pp — > fi fj, e + e~jj 
(or t t e + e~jj) and pp — > fi^e^bbjj (or pp — > r^e^bbjj) which break both lepton number 
and flavor and have no missing energy. The parent states could therefore be reconstructed 
up to the oringinal Wr decay. The lepton final states are predicted to have distinct flavor 
structures. We further point out that the l c — mixing also leads to the production of 
righthanded polarized top quarks from down-type squark decay, which is distinguishable from 
the X'QLd c trilinear couplings in the usual R-parity violating MSSM [14]. 

In section II, we study the general features in a class of models where R-parity is broken 
together with extra gauge symmetries. We derive new R-parity breaking terms from the 
Kahler potential and point out how to distinguish this class of model from others, e.g. the 
MSSM with usual R-parity breaking terms. In section III, we review the symmetry breaking in 
the context of minimal supersymmetric left-right (SUSYLR) model, emphasizing the necessity 
of spontaneous (dynamical) R-parity violation for SU(2)r gauge symmetry breaking. We 
discuss the flavor issues of R-parity breaking and its implications to neutrino mass and in 
section IV, we study the signatures of the model at the LHC. We mainly focus on the single 
production and decay of slepton NLSP via a heavier RH neutrino. The predicted multi-lepton 
final states and the flavor structure make the model distinguishable even in the early running 
of the LHC. Finally in section V, we point out a new contribution to the neutrinoless double 
beta decay in the model and conclude. 
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2 Spontaneous R-parity Violation with Extended Gauge Symmetry 



Unlike explicit R-parity violation, spontaneous R-parity violation (SRPV) has the advantage 
that it introduces only one new parameter into the R-parity conserving theory - the VEV 
of an R-parity-odd field. Furthermore, if R-parity violating scale is at the TeV range, above 
this temperature, R-parity is exact and therefore it is less constrained by cosmology. 

SRPV can be realized in various ways: in the first model where the idea was discussed [12], 
the superpartner of SM neutrino was given a non-zero VEV. Since lepton number is not a 
gauge symmetry of the MSSM, this leads to a doublet majoron which contributes to the 
.Z-boson width and LEP measurements therefore have ruled out this scenario. One could 
of course implement SRPV by the VEV of a right-handed sneutrino [13] in extensions of 
the MSSM that explain neutrino masses. Since the right-handed sneutrino field is a stan- 
dard model singlet, the majoron does not couple to the Z-boson and therefore escapes the 
constraints set by the LEP data. 

In this section, we will pursue the implications when the R-parity is spontaneously broken 
together with some extra gauge symmetry beyond Gsm = SU(2) lxU (l)y . Here we focus on 
the class of models where the extended gauge group G contains a subgroup SU(2)r. Clearly, 
the RH neutrino and its superpartner will be charged under the SU(2)r. Giving a non-zero 
VEV to u c will therefore give rise to new interactions. Models with dynamical R-parity 
breaking belong to this category. Furthermore, such model predicts that the scale of new 
gauge interactions is tied to the soft SUSY breaking scale. This leads to several interesting 
new features as we show below. 

The key distinguishing prediction of such a model is the existence of a large mixing 
between RH charged leptons and the gaugino superpartner of the Wr boson. Since we work 
with the gauge group G = SU(2)l X SU(2)r x U(1)b-l, the RH neutrino and charged lepton 
form a doublet under SU(2)r. 

After the RH sneutrino developing a VEV 
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it breaks both the SU(2)r gauge group as well as the R-parity, at the scale of (y c ). Due to 
the Higgs mechanism, the heavy gauge bosons acquire their mass by absorbing the scalars 
XmV c and l c as the longitudinal components. Due to supersymmetry, one would expect 
that a corresponding large Dirac mass would develop between the v c — Z'° and 1° — W^. 
Here we are interested in the chargino-lepton mixing, which is the new source of R-parity 
breaking effects. We explicitly write down the charged fermion mass matrix in the basis of 
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M 1/2 
M Wr rat 



(2.2) 
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where M\y R = 9r{v c ) is the Wr gauge boson mass and M\j2 ls the soft SUSY breaking 
mass for the chargino. The (1-2) element is absent because neither the RH neutrino nor the 
Higgs VEV couples l~ and W^. The determinant of this mass matrix is proportional to the 
(light) charged lepton mass m^ because the RH sneutrino VEV is an electroweak singlet and 
therefore does not break chirality. This means that there must be a physical state with the 
mass mi and identifiable as the charged lepton. 

Diagonalization of this mass matrix leads to a mixing between £ c and W^. This mixing 
Oiw is large if M x / 2 — M\y R . The physical charged lepton state is then given by 

t = 9 u i c + e m w+ , (2.3) 

where 0£w ~ 0(1)- We note there is no such mixing induced for W^. Due to this mixing, one 
can derive two new classes of R-parity violating interactions from the right-handed gaugino 
matter coupling terms: First from the gaugino coupling \/2gW^ v^i 1 ^ + h.c, we get 

Ci ew = V2g0 ew [iV e ,l' 4 + tvfl' "] . (2.4) 

Analogously, using the right-handed gaugino interaction with quarks and squarks, one can 
also write down the new RPV couplings for the squark-quark sector 

£9 ew = ^2 g0ew ^Vd c t + fctfd } . (2.5) 

Notice there is no supersymmetric counterpart of above terms generated, unlike those from 
the superpotential: XLLe c , X'QLd c , etc.. The new couplings we obtain here break not only 
R-parity but also supersymmetry, since we started from a mass matrix Eq. (2.2) including 
the SUSY breaking gaugino mass. 

As we will point out in Sec. 5, the most stringent constraints on the couplings in Eq. (2.5) 
are from neutrinoless double beta decay and HERA experiment [15], which tend to push the 
squark and gluino masses to TeV. On the other hand, the LEP2 Z-pole observables give a 
universal constraint [16] on the mixing parameter Qiw in both Eqs. (2.4) and (2.5), but it 
turns out to be rather mild. Therefore, the sleptons masses are still allowed to be not far 
above 100 GeV. 

From these new interactions derived from dynamical R-parity breaking, one would expect 
the following distinctive signatures at the LHC. 

• Single production of slepton NLSP via Eq. (2.4) and subsequent decays (see Fig. 1), 
which is the main topic being studied in this paper. 

• Top quark produced from down-type squark decay via Eq. (2.5), whose polarization is 
opposite to the components from the SM background, as well as the conventional A' 
term. 

We want to point out that these predictions are common to the models with R-parity bro- 
ken together with extended gauge symmetries, as well as to those breaking SU{2)r symmetry 
without the Higgs triplets, as long as the RH neutrino mass lies in the proper range. 
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Figure 1. Single production of RH slepton from RH sneutrino decays. The black box represents RPV 
e c — Wr mixing. 

3 An Explicit Model of Dynamical R-parity Breaking 

In this section, we present a model based on the group SU(2)l x SU(2)r x U(1)b-l where 
in the absence of R-parity breaking, the gauge symmetry does not break [2]. Thus the gauge 
dynamics dictates R-parity breaking. Hence it explicitly provides an example for dynamical 
R-parity breaking . 

The model considered in this section does not have the discrete parity symmetry. In the 
appendix, we will comment if such model can be built in a completely parity symmetric form. 
Here, we first review the salient features of the SU{2)l x SU(2)r x U(1)b-L model (without 
parity symmetry) for completeness and show how dynamical R-parity breaking occurs. 

3.1 SUSYLR: No gauge symmetry breaking without R-parity breaking 

The minimal SUSYLR model has the gauge group G LR = SU(3) C x SU(2) L x SU(2) R x 
U(1)b-l- The particle content and their representations under the gauge group for the 
completely parity symmetric case are listed in Table 1. The SU(2)r Higgs triplets A c , A c 
have been introduced to give mass to the RH neutrinos and facilitate the seesaw mechanism. 
For the sake of simplicity, we assume that the left-handed triplets A, A are decoupled at high 
scale and do not exist in the TeV theory. The superpotential of the model is 



where Y's are Yukawa couplings, / is the Majorana coupling of leptons and /xa is the //-term 
for triplets. Note that in the model there is no gauge singlets introduced. 

1 For recent papers where hidden dynamics breaks R-parity, see [17] and also some experimental implications 
of such models see [18]. However, in these works, the R-parity is not broken together with extended gauge 
symmetries that couple to SM fermions, and therefore, does not predict the phenomenology of dynamical 
R-parity breaking being discussed in this paper. 
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Table 1. Particle content in the minimal SUSYLR model. In this section, we first concentrate on the 
case with the left-handed Higgs triplets A, A do not exist in the TeV theory for simplicity. We will 
comment on the fully parity symmetric theory in the appendix. 



The corresponding soft terms are 



Koft = m~ (Q ] Q + Q C] Q C ) + mf (til 

+ X - (M 2L X a L X a L + M 2R X a R X a R + MxX bl Xbl + M 3 X g X g ) 

+ g T r 2 A^ i r 2 Q c + L T T 2 Alcf> t T 2 L c + iA f L cT r 2 A C L C 
+ BzatTr (T 2 4> T a T 2 4) b ) + B A Tr (A C A C ) + h.c. . 



+ L ct L c ) + m 2 A Tr(A ct A c ) + m\ Tr(A ct A c ) 



(3.2) 



The D-term potential as well as the scalar potential can be found in Refs. [2, 9]. 

The desired symmetry breaking pattern is SU(2) R x U(1)b~l — > U(l)y at the first step, 
giving definite meaning to the hypercharge Y = Is R + {B — L) /2, followed by the electroweak 
symmetry breaking. The key point to note is that the potential does not break any gauge 
symmetry in supersymmetric limit [19]. Even if the soft SUSY breaking terms are included, 
the gauge symmetry still remains unbroken as long as the RH sneutrino has zero VEV. Parity 
and SU(2) R x U{1)b-l — > U(l)y breaking become possible only if the RH sneutrino picks up 
a non-zero VEV. The RH sneutrino being superpartner field has odd R-parity and therefore 
its VEV breaks R-parity - hence the claim [2] that there is no parity breaking without R- 
parity breaking in the minimal SUSYLR model. Furthermore, it was shown in [6] that RH 
sneutrino VEV is tied to the soft mass scale Mstjsy • This implies an upper bound on the Wr 
gauge boson mass of order of the SUSY breaking scale. 

To see this explicitly, we write down the potential including all VEV's given below; 
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(3.3) 



where we choose to break the R-parity along the RH electron sneutrino v% direction, for 
phenomenological consideration to be explained in Section 3.2. The scalar potential involving 
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Figure 2. Correlations among the VEVs. The left-panel tells us that R-parity is broken as much as 
parity, (y^) 2 ~ vj^ + vj^. The middle panel shows that the minimum always points towards the flat 
D-term potential direction, (vl) 2 = 2(v 2 i — vj i ). The right panel tells the value of the VEV is related to 
parameters in the potential as (yf) w A/if 2 . This agrees with the upper bound obtained in Ref. [9]. 



the Higgs triplets and RH sneutrinos is 

V = M 2 v 2 R + M 2 v\ - 2Bvrvr + \f\ 2 (u c e ) 4 + [4|/| 2 4 + m 2 - 2\A\v R - 2\f\fi A v R ] {V c e ) 2 

+ \{g 2 +g' 2 )m) 2 -2vl + 2v R ? , (3.4) 

where M 2 = /i 2 ^ + m 2 ^ , M| = /i 2 ^ + m 2 ^ and B = B/± . For simplicity, we have assumed the 
matrices /, Af, m~ to be flavor diagonal and / = f ee , A = (Af) ee and 771q = (m~) ee . The 
VEVs of the Higgs bidoublets have been neglected in the first stage of symmetry breaking. 
The potential should satisfy B < M1M2 to be bounded from below. This can be seen by 
considering D-flat directions (z/g) = 0, (A) = (A c ) = vt\ and (A) = (A c ) = vti, where t\ is 
the Pauli matrix. This scalar potential has the property that on (i/g) = surface, there is no 
symmetry breaking, i.e., v R = v R = at the minimum. The acceptable minimum that breaks 
parity and the gauge symmetries therefore necessarily breaks R-parity. 

We have carried out numerical study of the minimization of the scaler potential, by 
scanning over the bulk of parameter space (Mi, M2, vB, A, /^a> mo) £ [100, 1000] GeV and 
/ G [0.1, 0.5]. The true vacuum must satisfy V m i n < and (uf) 7^ 0. It turns out that there are 
interesting correlations among the VEVs of RH neutrino and the Higgs fields. They are shown 
in Fig. 2. Typically, we find the D-term potential always vanishes, i.e., (z?g) 2 = 2(v R — v R ). 
Therefore, the physics at the RH scale does not bring additional terms to the Higgs potential. 
The sneutrino VEV and the Higgs triplets VEVs are of the same order, (y^) 2 ~ v R + v R , as 
well as an approximate relation (u^) ~ A/3/ 2 . This agrees with the upper bound obtained 
in Ref. [9]. The key point is that, in SUSYLR model, the right-handed scale is dynamically 
generated through the SUSY breaking soft mass scale [2] , 

vr < p , (3.5) 
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where -Msusy ~ 0(100) GeV corresponds to the generic soft SUSY breaking mass scale. 



Vff /, VR, m (GeV 4 ) with <;/> nii „=20S5 GeV 




v/?(GeV) v R (GeV) 

Figure 3. Left panel: The potential V as a function of vr for given (yf) — 2055 GeV, vr — 1607 GeV 
(solid curve) and (yf) = vr = (dashed curve). Right panel: Contour plot of the potential V in the 
vr — vr plane for (u^) = 2055 GeV. 

In order to illustrate the role of R-parity violation in symmetry breaking, we choose the 
following set of parameters 

Mi = 213 GeV, M 2 = 251 GeV, VB = 150 GeV, 
/iA = 517 GeV, A = 240 GeV, m = 376 GeV, / = 0.21 . (3.6) 

The resulting VEV's and the minimum potential value are 

(v c e ) = 2055 GeV, v R = 2063 GeV, v R = 1607 GeV, V mill = -1.0 x 10 12 GeV 4 . (3.7) 

The configuration of the potential around the vacuum is shown in Fig. 3. Clearly, the global 
minimum of the potential breaks R-parity, i.e., (u^) ^ 0. Because the R-parity and the 
lepton number are broken simultaneously with the gauge symmetries, no massless Majoron is 
present. On the other hand, the dynamical R-parity breaking associated with gauge symmetry 
breaking at few TeV scale offers rich phenomenology. 

3.2 Neutrino mass and flavor alignment of R-parity violation 

Neutrino masses in this model have been discussed extensively in [9]. We review the salient 
points for completeness and, in particular, constraints on flavor of the R-parity violation. In 
SUSYLR model, the matter fields obtain their Dirac masses from the coupling to the Higgs 
bidoublets $1,2- Generally, there are four SU(2)l Higgs doublets at the electroweak scale. 
The additional neutral Higgs bosons will lead to flavor changing neutral currents at tree level. 
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This can be suppressed either by proper doublet-doublet splitting or by cancellations [20]. 
Another way to avoid the flavor changing Higgs effects is by replacing the second B — L = 
bidoublet with a B — L = 2 bidoublet as has recently been suggested [21]. We do not discuss 
this here. 

We assume that the bidoublet Higgs fields take the following form of VEV's, 









, $2 = 


K2 

























(3.8) 



where K\ gives Dirac masses to the up-type quarks and neutrinos, while k 2 contributes to 
down-type quarks and charged leptons masses. With this VEV structure, the Wl — Wr 
gauge bosons do not mix with each other. Here we mainly focus on the lepton sector. We 
will attribute the hierarchies among the charged lepton masses to the Yukawa couplings. In 
particular, we focus on the low tan/3 = k 2 /ki ~ 0(1) regime. In this case, the Yukawa 
coupling constants are set as y T « 10 -2 , y M « 10 -3 , y e « 10 5 and {Yv)ij ~ 10 -6 . Even 
though there are four SM Higgs doublets (or two MSSM Higgs pairs), since only two of them 
contribute to fermion masses and the other two play the role of spectators, our tan/3 is same 
as the MSSM one. The /i$ and B§ play a similar role as the \i and B^ parameters in the 
MSSM for the electroweak symmetry breaking. 

Because of R-parity violation, there are additional contributions to neutrino masses, on 
top of the Type-I seesaw mechanism. They arise from neutrino-neutralino mixing which has 
been calculated in Ref. [9] with R-parity breaking in the RH electron sneutrino direction, 



(M2L M2R Ml\ 

2M 2L M s ( f i 11 ti 22 - tf 2 ) V 9l Sr Sbl/ 



(3.9) 
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^> 2 



assuming the LH sneutrino VEV's are negligible. Ml,r and Mbl are soft supersymmetry 
breaking gaugino masses. The corresponding Feynman diagram is shown in the left panel 
of Fig. 4. Choosing the R-parity violation along the v\ direction helps to avoid too large 
contributions through the y T , y M couplings, while the electron and neutrino Yukawa couplings 
are sufficiently small to keep the neutrino mass scale 0(0.1) eV in tact. 

We also notice that there are radiative corrections to the neutrino mass [22, 23], which 
are also proportional to the corresponding Yukawa coupling and are loop suppressed (see the 
right panel of Fig. 4, as well as Fig. 2 in [24]). So they are safely small as long as (ufy and 

are vanishing and tan j3 is low. 

The above discussions justify our choices of VEV configuration in Eq. (3.3). 



Before closing this section, we comment on the flavor violations. The next section will 
mainly concentrate on the scenario with slepton NLSP singly produced from the gauge 
boson resonance, hence we need to understand the existing experimental constraints on the 
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relevant mass scales. In the minimal non-SUSY left-right model, the famous neutral K- 
meson mixing tends to push to Wr mass to be above 2.4 TeV [25]. In the supersymmetric 
version, loop diagrams mediated by superpartners also make additional contribution to both 
quark and lepton flavor violation processes. They are safely small if the relevant mass scales 
are high enough. Otherwise, in order to optimize the discovery prospects at the LHC, the 
superpartners and W/j-boson masses have to lie in the (sub-) TeV regime, which requires fine- 
tuning the flavor structures of the model in a similar way to the MSSM. This is nothing 
but the SUSY flavor problem, and the constraints on scales are quite model dependent. 
In principle, there could also be contribution to flavor violations from higher dimensional 
operators controlled by unknown physics in the UV. 

Therefore, in the following we shall only adopt the bounds M\y R > 1 TeV and mj > 
100 GeV from Tevatron and LEP2 direct searches, respectively. 

<v°> <v°> 
X X 

I I 
y 1 gK m gK t y 

— > 1 a • s >< " ' a 1 « — 
v h v v h v 

Figure 4. Contributions to neutrino masses from R-parity violation. Left panel: tree-level contri- 
bution due to neutrino-neutralino mixing. Right panel: loop-suppressed radiative correction to the 
neutrino masses. The v represents neutral gaugino fields. The black dots are the usual Higgs VEV 
insertions. 
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4 Single Production of Slepton NLSP and its Decay at the LHC 

In this section, we start exploring the LHC implications of this model. First we need to know 
the R-parity violating (RPV) interactions that induce the decays of the sparticles produced. 

4.1 Relevant RPV couplings 

In general, spontaneous R-parity breaking through the RH sneutrino VEV generates the 
bilinear terms in the superpotential and the soft potnetial 

Wjf= fii.LiH u , V>„ soft = BiLiH u + h.c. . (4.1) 

The bilinear term facilitates the R-parity breaking decay of the lightest neutralino, xj as 
follows: x? ~~ * Z°v, Xi W^lzp or x5 — > ^\^2 V [^]- I n the literature, more complete collider 
phenomenologies of R-parity violation from the superpotential has been studied in detail and 
reviewed in Refs. [22, 26-28]. 
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In the SUSYLR model, the bilinears arise from the electron and neutrino Yukawa cou- 
plings and the corresponding A-term once the RH sneutrino VEV is inserted [9]. Therefore 

— - - y e , - i(T 6 . (4.2) 

These bilinear terms will induce trilinear R-parity breaking terms XLLe c , X'QLd c . The most 
important terms for the following study are those associated with third generation fermions, 

X\Uib c + X l ul l r c , (4.3) 

where A^ = ytfJ-i/n®, X% = Vr^ij V* an d i = 1,2, 3. On the other hand, the A" term will not 
be generated, since baryon number symmetry is respected by the sneutrino VEV, thereby 
guarantee the proton stability. 

As already stated in previous sections, a distinct feature that arises when R-parity is 
dynamically broken together with an SU{2)r gauge symmetry is the large mixing between 
the RH electron e c and the SU{2)r gaugino, i.e., Wj£. This is not present in the MSSM with 
general R-parity violation. The form of charged fermion mass matrix and the obtained mixing 
in SUSYLR model is given explicitly in the Appendix. From the usual gaugino Yukawa-like 
coupling term for fx and r flavors, one obtains the new couplings (similar to Eq. (2.4)) 

4 ew = V2g9 eW [e c ^p + e c v°Ji c ] + V2g6 eW [e c i^f c t + e c ^r c ] . (4.4) 

As we will see, these interactions open new channel for the single production of a slepton at 
hadron colliders (Fig. 1). Similarly, from the neutralino mass matrix, one can also obtain 
large mixing between the Z' and u°, and in turn the couplings 

££ w = V2g z , 6 NZ , Uif~^ + + + v^f] + h.c. (4.5) 



In principle, sparticle single production could also happen through the mixing between 
and Z' [24] which, however, calls for some tuning between Mz> and M\/2- 

Contrary to the usual R-parity breaking term from superpotential, these new R-parity 
breaking sources come from the gaugino Yukawa- like couplings (in the Kahler potential). As 
we illustrate in the below, such theories could be tested at the LHC where the new gauge 
interactions are accessible. 



4.2 Branching ratios of slepton NLSP decay 

From the previous sections one learns that, in the SUSYLR model under discussion, a new 
class of RPV couplings Eq. (4.4) emerge due to the mixing between e c and W^. To study 
its implications for hadron colliders, we need to know the sparticle spectrum. If one takes 
the assumption of universal scalar masses at high scale, the RH sleptons are likely to be the 
lightest among matter superpartners in the MSSM due to the smaller Yukawa couplings as well 
as the smaller weak gauge couplings [29] . The situation would be similar in SUSYLR models. 
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The lightest slepton could be a stau or the smuon depending on detailed parameter range. 
In our study, we will assume that smuon is the lightest superpartner above the gravitino, the 
latter in our model could be the very weakly unstable dark matter. 

As promised, we study the implications of scenario at the LHC for the case where smuon 
or stau is the NLSP among the superpartners. Due to the relatively low tagging efficiency of 
the tau lepton, we would focus on the smuon. 

The new LHC signals originate from the production of Wr in pp collision and its subse- 
quent decay to muon and RH muon neutrino which subsequently decays. In the non-SUSY 
LR models with type I seesaw the RH neutrino decays mostly to the three body final state 
i^i^jj via Wr exchange [30, 31]. However in the SUSY version, if the smuon, Jl is lighter 
than v c , an interesting new two body final state channel emerges: RH neutrino decays to a Jl 
and an electron. Since this is two-body decay, for the smuon sufficiently light, it will certainly 
dominate over the three body non-SUSY mode, 2 as shown in Fig. 5. Therefore, the smuon 
single production could take this advantage and be large enough to be probed at the LHC. 




Figure 5. Production of the RH neutrino and its decays. Left-panel: the usual same-sign lepton 
diagram for Wr discovery. Right-panel: Single production of the smuon through RH neutrino RPV 
decay. In this case, there is equal possibility to break the muon lepton number twice through the 
Majorana mass of v° or not, so one can get either Jl + e] i or Jl~ from its decay. The black box 
represents sneutrino VEV insertion as indicated in Fig. 1. 

As the NLSP, the smuon /J NLSP = cos a JIr + sincve^jui, will decay dominantly through 
R-parity breaking interactions rather than the Planck scale suppressed decay to the gravitino. 
In the case where there is a large mixing between LH and RH smuons, i.e., sin a ~ 0(1), Jl + 
can decay to tb or tv through the induced trilinear RPV terms as shown in Eq. (4.3), although 
suppressed by the small y e or y v Yukawa couplings. Note that electroweak symmetry forbids 
the direct coupling of RH smuon to //, and the RH sneutrino VEV does not help because it 
is also a singlet. Only the LH and RH smuon mixing term (rn~)LR which is proportional to 
the Higgs doublet VEV, can facilitate this decay. 

On the other hand, if the mixing term (m~)LR is severely suppressed, i.e., sin a <C 1, the 
smuon is almost purely RH. Therefore, it decays through a four-body channel, as shown in 
the right panel of Fig. 6. Such decay rate is proportional to the gauge coupling instead of the 

2 Since the RH neutrino mass matrix is proportional to the matrix /, which we have taken to be diagonal 
in the basis of physical charged leptons, there is no further flavor changing in v a mass matrix (propagator). 
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Figure 6. Two- and four-body decay modes of the smuon NLSP. The black box represents sneutrino 
VEV insertion as indicated in Fig. 1. The black dot stands for the usual Higgs VEV insertions. 
Hereafter, we denote the RH smuon as Jl c = J2~^. 

small y e or y u couplings. It could be comparable or even dominate over the above two-body 
decays when the latter is further suppressed by the LR smuon mixing. In principle, Jl c can 
decay to both e + u^jj final states. However, since the intermediate RH neutrino is off-shell, 
the probablity to break the muon lepton number is larger than that conserving the lepton 
number. This point can be see from the blue and brown curves in the left panel of Fig. 7. 

The branching ratios to different final states of the smuon decay has been plotted in Fig. 7 
with the following parameters chosen, M\y R = 2TeV, M u c = 500 GeV and g6 e w = 0.2. In 
the suppressed LR slepton mixing case (left panel, here and below, we will take sin a ~ 10 -3 
as a benchmark point), the four-body decay of smuon NLSP always dominates over all the 
two-body channels. In the large mixing case (right panel), the tb channel will dominate if 
it is kinematically allowed, while tU and the four-body channels e + fi~jj could respectively 
dominate in certain low smuon mass windows. In both cases, since the Majorana RH neutrinos 
are involved in the production and/or decay processes, lepton number can be broken, which 
leads to the most promising discovery channels at the LHC. The expected signatures are 
listed in Table 2. 





large Jl — Jl c mixing 


suppressed mixing 


M^ c > m t + mi, 


pp — > n^e^tb, fi+e + tb 


pp — > (j, fj, e + e~jj 


Mjic < mt + mb 


pp — > [i^fi e + e~jj 
pp — > fi~e + T + + J^t 



Table 2. Expected final states in the single production of the NLSP /i c with M^c < M„c assumed. 
Large and suppressed Jl — Jl c mixing cases are both listed. 

In the following two subsections we will discuss the signature of the single production and 
decay of smuon NLSP via a heavier right-handed neutrino. We also discuss possible standard 
model backgrounds and elaborate on the selection criteria necessary for such signals to be 
significantly observed over the standard model background. The large number of diagrams 
involved in the standard model background processes are calculated using the helicity am- 
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plitude package MadGraph [32] and CalcHEP 2.5.4 [33]. To estimate the number of signal 
and background events as well as their phase space distribution(s), we use a parton-level 
Monte-Carlo event generator. In our numerical analysis, we use the CTEQ6L parton distri- 
bution function [34] and fix the factorization scale Q 2 = s/4. In our parton-level simulation 
of both signal and background events, we smear the leptons and jet energies with a Gaussian 
distribution according to 

AK a 

b (4.6) 



with the CMS parameterization, ai = 5%,bi = 0.55% and aj = 100%, bj 
sum in quadrature. 



Vo, © denotes a 



4.3 pp 



e+e A i± /x ± jj 



This particular final state dominates when the mixing between LH and RH smuon is sup- 
pressed (or in the low mass (~ M top ) region for a large mixing). In this section, we will denote 
smuon NLSP as Jl c since it is mainly the RH component. The most striking feature of this 
final state is the three same sign leptons and one opposite sign lepton associated with two 
jets without missing energy. Assuming the narrow width approximation for u c and £t c , we 
can simply write down the signal cross-section o~ s (pp —> e + e~ fi + j j) as 



a(pp — > e + e /U + /i + jj) « a(pp — > W. 
x Br(f£->/2 c e~) x Br(/2' 



r -> 

-> e Vjj) + Br(i£ 



(4.7) 



jj d *e + ) x Br(/i cT ^e fi + jj 
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Figure 7. Branching ratios for the smuon NLSP decay. The left panel represents the suppressed LH 
and RH slepton mixings (~ 10 -3 ) case, while in the right panel, we take an unsuppressed 0(1) such 
mixing. Charge conjugated final states are not listed but also possible. 
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where the red (dotted) arrow indicates lepton number violation by two units on the involved 
RH neutrino propagator. The charge conjugated final state cr(pp — > e + e~ /U~ pT j j) which is 
mediated by the intermediate boson can be similarly approximated. In our analysis, we 
combine both these two final states. We define the signal identification with four charged 
leptons and two jets. The events are further selected by the following set of cuts 

1. We require that both jets and leptons should appear within the detector's rapidity 
coverage, namely 

\ V (£)\ < 2.5, \ V (j)\ < 3 . (4.8) 

2. The leptons are ordered according to their transverse momentum (pr) hardness and the 
Pt of the leading lepton must satisfy 

PT (h) > 100 GeV , (4.9) 

and for rest of the leptons 

p T (£) > 15 GeV . (4.10) 
For two associated jets we demand that 

f^ ts >25GeV. (4.11) 

3. We must also ensure that the jets and leptons are well separated so that they can be 
identified as individual entities. To this end, we use the well-known cone algorithm 
defined in terms of a cone angle AR a p = J {Acj) a p) 2 + {Ar] a p) 2 with Ac/) and Arj being 
the azimuthal angular separation and rapidity difference between two particles. We 
demand that 

AR j:j > 0.4 , AR ej > 0.4 , AR U > 0.4 . (4.12) 

4. In our analysis, We use simplified definition for the missing transverse energy: I$ T = 

\j i^Px) 2 + (52PyY > where the sum goes over all observed charged leptons and jets. 
We demand that there is no significant missing energy in our signal 

Tfij. < 30 GeV. (4.13) 

Our choice of px cut on the leading lepton (Eq. (4.9)) can be well justified from the pt 
distribution of all four leptons as displayed in Fig. 8 assuming M\y R = 1 TeV, M v c = 500 GeV 
and M^c = 300 GeV and at -y/s = 14 TeV. Here, one should note that while generating pt 
distributions (Fig. 8), we impose an uniform loose cut (pt > 15 GeV) on all four leptons, 
however, rest of the cuts remain unchanged. From the choice of mass parameters and simple 
kinematics of the production and decay chain, it is very obvious that the leading lepton (£i) 
comes from the two body decay of heavy W% — > fi + + v c , while rest of the leptons originat- 
ing from the cascade decay chain of u c and jl c have relatively softer transverse momentum 
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Figure 8. pt distributions of all four leptons in the process pp — > e + e _ /Lt ± /x ± > jj at ^fs = 14 TeV. The 
leptons are ordered according to their px hardness (pr{h) > Pt(^) > Pt(^) > Pt(^4))- We have 
fixed M Wr = 1 TeV, M u , = 500 GeV, and = 300 GcV. 
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Figure 9. Signal cross sections a(pp — > e + e~ fJ^^jj) (after all cuts as mentioned in the text) as a 
function of smuon mass at the LHC with ^fs = 7 TeV, 10 TeV and 14 TeV. Three curves from top 
to bottom in each panel correspond to Mw R — 1 TeV, 1.5 TeV and 2 TeV respectively. M u <= is kept 
fixed at 500 GeV. 

compared to the pt of the leading lepton. On the other hand, as the RH neutrino mass is 
increased to a value closer to the Wr mass, the first lepton becomes softer. However, in this 
case, the lepton from the decay — > e~Jl c merits the highest px and will serve as the hardest 
lepton (^i). 

In Fig. 9 we show the total signal cross-section a s (after imposing all the cuts mentioned 
above) for the process shown in Eq. (4.8), as a function of the smuon Jl c mass at the LHC for 
7 TeV, 10 TeV and 14 TeV energies. In each panel, three curves from top to bottom correspond 
to My/ R = 1 TeV, 1.5 TeV and 2 TeV respectively. We fix the RH neutrino mass M u c = 500 
GeV and the mixing parameter gO e w = 0.2 for the present analysis. Before estimating the 
possible Standard Model backgrounds to this particular channel, we would like to discuss the 
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general behaviour of the signal cross sections. 



• In all three panels, irrespective of Myy R , the a s first rises with the increases of smuon 
mass and then becomes almost flat and finally drops sharply as M-^c becomes degenerate 
with right-handed neutrino mass M v c. 

• The initial rise of the cross-section with the smuon mass can be understood from the 
fact that for lighter smuon mass (M^c ~ 100 — 200 GeV), the decay products of smuons 
Jl c — > e + fi + jj are more collimated and fail to satisfy our isolation criteria for leptons 
and jets as shown in Eq. (4.12). As the smuon mass increases, leptons and jets which 
originate from the cascade decay of smuon tend to appear with larger AR, thus satisfying 
the isolation criteria as displayed in Eq. (4.12). As a consequence, the a s for heavier 
smuon mass (M^c < M v c) is significantly larger than for lower smoun mass region. 

• The signal cross secion a s strongly depends on ^/s, mass Mw R and off course on M^c . 
There is a possibility that the LHC may also run at \fs = lOTeV, before attaining to 
its designed y/s = 14 TeV. Keeping this in mind, we decided to provide our observation 
for y/s = 10 TeV also. In is very interesting to note that for all the choices of M\y r and 
\/s the smallest cross-section always correspond to M^c = 100 GeV, while the largest 
one correspond to M-^c which lies between 400 — 430 GeV as shown in Table 3. 



M Wr 


7 TeV 


10 TeV 


14 TeV 


ITeV 


3.2-46.8 


7.0-100 


13-178 


1.5 TeV 


0.3-4.7 


0.9-14.6 


1.0-31.7 


2 TeV 


0.035-0.5 


0.1-2.3 


0.4-6.2 



Table 3. The range of minimum and maximum cr(pp — > e + e~/j, fj, jj)(fb) at the LHC for y/s — 
7, 10, 14TeVandMvv H = 1, 1.5, 2 TeV, respectively. The corresponding smuon masses are mentioned 
in the text. The other parameters are taken as M v o = 500 GeV and g0 e w = 0.2. 

Mass reconstruction: The most important feature of our signal events is the effective 
reconstruction of all three heavy particle masses from the final state charged leptons and 
jets. We first select two softest leptons (satisfying our selection criteria) from the four lepton 
set and then recombine these two leptons with the two jets to reconstruct the smuon mass, 
Mjj£ 3 i 4 ~ M^c After the obtaining the smuon resonance, we attempt to reconstruct the RH 
neutrino mass by combining two jets, two softest leptons with one of the two hardest leptons 
l\ or 1%. In this case, we face the complication due to combinatorics with two choices of 
pairing for with Mjj£ 3 £ 4 . Finally, Wr can be reconstructed by combining all four charged 
leptons and two jets. In Fig. 10, we display the invariant mass distribution for Jl c , and Wr 
at 7 TeV LHC. Fitting the mass distribution with a Gaussian, we get the following values 

MS = 301.36 ± 1.74 GeV, M$ = 500.43 ± 0.75 GeV, M$ R = 999.45 ± 2.10 GeV ,(4.14) 
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Figure 10. Invariant mass distributions for M^c, Af^c and Mw R in the pp — > e + e fi (Jt jj process at 
= 7 TcV with 3 ftT 1 data for Afy = 300 GeV, M^ r = 500 GeV and M Wr = 1 TcV respectively. 
The error-bars shown are statistical only for the indicated luminosity. Results of Gaussian fitting are 
also shown. 

where the input masses considered for this mass reconstruction procedure are the following 
M| c ue = 300 GeV, M*f e = 500 GeV, M^ e = 1000 GeV . (4.15) 



SM background 


cr (pb) 


a e ± e ± (fb) 


pp — > bbbb 


387.5 


0.16 


pp — > ti 


448 


0.09 


pp — > Z°bb 


0.051 


3 x 10" 5 


pp -> W^^^Z 


6.7 x 10" 4 


2 x 10" 5 






0.25 



Table 4. The list of leading-order SM backgrounds that could mimic our signal. <jq and <j£±i± are 
defined in the text. These numbers correspond to y/s = 14TeV. 

SM backgrounds: In principle, there is no intrinsic standard model background to the 
AL = 2 processes. However, there some standard model processes which could mimic our 
signal if the missing transverse momentum of neutrinos are balanced. One of the dominant 
background is pp — > bbbb, followed by semileptonic decay of all the b-quarks. We generate this 
background using with the following basic cuts pr(b) > 25 GeV, \r](b)\ < 2.5 and AR^b > 0.4. 
The leading order cross-section is 388 pb at ^fs = 14 TeV. After hadronization, one of the 
B° or B has to oscillate before decay, in order to get a pair of same-sign dileptons. The 
probablity of having bb — > e^/i 1 * 1 , / u =t ^ =t is about P^±^± ~ 2 x 10~ 5 , as estimated in [35]. After 
taking into account the semileptonic branching ratio ~ 10% for the other two 6-quarks we find 
this background cross-section ~ 10 _1 (fb). The other aparently looking very severe standard 
model background is pp — > ti. At leading order, the top pair production cross-section (a t i) 
is 448 pb at the LHC with ^/s = 14 TeV. After taking into account the leptonic branching 
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fraction of two W bosons (from t — > bW + ) and P»±e±, the rate goes down to ~ lO^ 1 (fb). 
Here, we would like to mention that if we take into account the higher order QCD effects, 
a t i becomes ~ 900 (pb) , which means our final background cross-section from ft process may 
increase atmost by a factor of two. The other sub-leading standard model processes which 
my fake our signal processes are pp — > Zbb, pp — > W^W^W^ Z° and pp — > W^W^W^h. 
In the case of pp — > Zbb, process, Z — > £ + £~,£ = e,[i and same sign leptons will come 
from bb pair by oscillation of one of the B° meson before decay. As a result of this, the 
a(pp —> Zbb) will be suppressed by Bi(Z — > l + £~),£ = e,fi and Pj?±e±- The rate for same 
sign leptons from remaining two processes are negligibly small. In Table 4, we summarize the 
standard model background cross-sections, where, oq and o-£±i± correspond to the leading 
order cross-sections before and after folding with different suppression factors arising from 
leptonic branching ratios of W^,Z bosons, semi-leptonic branching ratio of b(b) quark and 
finally Pj!±£± respectively. From this very simple minded exercise, we conclude that our signal 
is almost SM background free. 



4.4 pp — > ^ ± e bbjj 

The smuon heavier than top quark and with large mixing L-R mixing 0(1) can lead to this 
final state. Here we call the smuon NLSP as ju, without definite chirality. As shown in Table 2, 
the smuon will dominately decay to tb via the A' coupling in this case. For this signal topology, 
we select events with two same sign different flavoured (SSDF) charged leptons and four jets. 
The cross section for this channel in the narrow width approximation can be expressed as 

a(pp — > n + e + bbjj) « a(pp —t W£ — > ^ + v c ^) ■ Br(t^ — > fl~e + ) ■ Br(/i~ ->■ tb) ■ Bi(t — > bjj) . 

(4.16) 

The signal also includes the charge conjugated final state cr(pp — > fi~e~bbjj) via intermediate 
boson. 

Our selection cuts are same as shown in Eqs. (4.8) — (4.13), except for the transverse 
momentum cut on the jets. After ordering all four jets according to their p?, we impose 
following cut on the hardest jet (ji): 

PT (ji) > 60 GeV (4.17) 

and for rest of the jets 

Pr(i2,i3,i4)>25GeV. (4.18) 

In Fig. 11 we display the px distribution of two leptons and four jets respectively after 
ordering them according to their px- While generating these distributions, we impose the 
following cuts on the pt of leptons and jets, rest of the cuts remain unchanged, 

p T (^)>10GeV, p£ ts >15GeV. (4.19) 
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P, (GeV) Pt (GeV) 

Figure 11. pt distributions of two leptons (left-panel) and all four jets (right-panel) in the process 
pp — > ^ ± e ± bbjj at ^/s — 14TeV. The leptons and jets are ordered according to their px hardness 
{Pt{Zi) > Pt{£i) and (pr(ji) > Pt{32) > PtHz) > Priji)- The other model parameters are same as 
in Fig.8. 

We take the same set of mass parameters as in the previous subsection. In this case too, the 
leading lepton comes from the two body decay Wr —> ii + v^. On the other hand, the leading 
jet ji mainly comes from the two body decay of the smuon, while the second hardest jet is 
produced from the top quark decay. From the nature of the pt spectrum of leptons and jets 
as shown in Fig. 11, we can justify our choice of px cuts (Eqs. (4.17) and (4.18)) used in this 
analysis. 

In Fig. 12 we show the signal cross section (after all cuts on final state leptons and jets 
as mentioned above) for this channel as a function of the smuon mass at the LHC for 7 TeV, 
lOTeV and 14 TeV energies. In each panel, three curves from top to bottom correspond to 
M Wr = 1 TeV, 1.5 TeV and 2 TeV respectively. M v o is kept fixed at 500 GeV and g6 eW = 0.2, 
the same as in Fig. 9. Comments on the cross sections are in order. 

• In this case, since we look for Jl + — > tb, we focus on the smuon mass above the top quark 
threshold, as is displayed in all three panels of Fig. 12. As the smuon mass increases, 
the leptons and jets originating from the cascade decay of smuon tend to appear with 
larger AR between each other, satisfying the isolation criteria shown in Eq. (4.12). 

• The signal cross section begins to drop for heavier smuon mass (> 350 GeV) irrespective 
of M\\r R and choice of the LHC energy. This is mainly due to the branching ratio 
suppression of the Jl + — > tb decay mode, as can be seen in the right panel of Fig. 7. 
Secondly, there is also the phase space suppression when becomes close to right- 
handed neutrino mass M v c 

• In Table 5, we show the range of signal cross sections for different values of at the 
LHC with = 7, 10, 14 TeV and M Wr = 1, 1.5, 2 TeV, respectively. The other 
parameters are taken as M u c = 500 GeV and gQ e w = 0.2. We quote the minimum and 
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Figure 12. Signal cross sections cr(pp — > fi~e~bbjj) (after all cuts as mentioned in the text) as a 
function of smuon mass at the LHC with ^fs = 7 TeV, 10 TeV and 14 TeV. Three curves from top to 
bottom in each panel correspond to Mw R = 1 TeV, 1.5 TeV and 2 TeV respectively. M v c is kept fixed 
at 500 GeV. 



maximum values of the signal rate. For all three choices of M\y R and yfs the smallest 
cross section always correspond to the value of M-n which is colse to M u c : while the 

r* j-L 

largest cross section correspond to lying between 260 — 300 GeV. The signal gets 
enhanced by more than factor of 2 as the LHC energy increases from 7 TeV to 10 TeV, 
and by another factor of 2-3 up to 14 TeV. 

Mass reconstruction: We now discuss the mass reconstruction strategy of all three heavy 
particles from the final state charged leptons and jets. From the sample of four jets, the 
hadronically decaying SM TV-boson is reconstructed from pair jets whose invariant mass 
{rrijj) is closest to Mw- The top quark is then reconstructed from the reconstructed W and 
one of the two remaining jets. We select the one which gives a invariant mass closest to Mt. 
The smuon mass is reconstructed from this Mt and with the last jet = my. Next, we 
attempt to reconstruct the right-handed neutrino mass by combining with one of the two 
leptons i\ or l<i- In this case, we are facing the combinatorical background with two choices 
m tje 1 , mijf 2 - Finally, the W/j-boson mass can be reconstructed by combining all four jets and 
two charged leptons. We will not explicitly show the reconstruction figure here, which looks 
very similar to Fig. 10. 

SM background: In this case, the standard model process which can mimic our signal is 

pp tiW ± bbW + W~W ± jjbb^i^ , (4.20) 

where £,£' = e,/U. In our analysis, we do not impose the requirement of b tagging, since the 
standard model background also contains 6-jets, and b tagging would not improve the signal 
significance considerably. The standard model background cross sections from pp — > tiW^ 
process is shonw in Table. 6, at different LHC energies. We expect this rate would further go 
down significantly (by several orders of magnitude) once we impose our selection criteria on 
the final state leptons and jets. 
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M Wr 


7 TeV 


10 TeV 


14 TeV 


ITeV 


0.3-21.22 


0.78-47.9 


1.4-85 


1.5 TeV 


0.18-3.04 


0.58-9.6 


1.3-21.32 


2TeV 


0.05-0.38 


0.22-1.75 


0.62-4.73 



Table 5. The pp — > fi ± e ± bbjj signal cross sections (in fb) at the LHC for y/s = 7, 10, 14 TeV 
and M\y r = 1, 1-5, 2 TeV, respectively. The other parameters are taken as M„c = 500 GeV and 
g^eW = 0.2. Here we quote the minimum and maximum values of the signal rate and the corresponding 
smuon masses are shown in the text. 





u tiw ± (fb) 




7 TeV 


99 


1.32 


10 TeV 


206 


2.77 


14 TeV 


377 


5.05 



Table 6. The dominant SM background pp — > ttW ± that could mimic our signal. <T t tw ± correspond 
to the production cross-section of ttW ± and represents cross-section for bbjj/d e final state 

before any cuts. 

We also comment on the other standard model background pp — > bbjj, which has a huge 
cross section ~ 10 5 pb after basic cuts. Taking into account of the oscillation of bb to get 
same-sign e^/i 1 * 1 Pe±e± wu l reduce it down to the order of ~1 pb. The cuts on missing energy 
and the hardest lepton and jet will further reduce the cross section. Moreover, in this case, 
highly energetic 6-jet will produce charged leptons which will be very close to the associated 
c-jet, as a result of this, lepton-jet isolation criteria will play a decisive role in reducing this 
background further. Therefore, we conclude that this background will be also under control. 
The remaining backgrounds pp — > W^W^W^ Z° , pp — >• W^W^W^h and pp — > jjjjW^W^ 
are much smaller [36]. 

5 Some Generic Low-energy Constraints 
5.1 Neutrinoless double beta decay 

In this model, there are several new contributions to neutrinoless double beta decay in addition 
to the usual light neutrino contribution. The contribution from the RH neutrino exchange as 
in the non-SUSYLR models was already discussed [37] . 

In our model, there are two new contributions arising from the e c — Wr mixing. The first 
one is given in left panel of Fig. 13 below. Its contribution to the effective neutrino mass is 
given by 

\M Wr J m w 
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where pf ~ 50 — 100 MeV is the typical momentum transfer in this process. For 8 e w ~ 0(1), 
it is of same order as the RH neutrino contributions to this process in non-supersymmetric 
case. 

The second contribution is given in the right panel of Fig. 13 and the corresponding 
effective neutrino mass is 



Note that for this to be consistent with the current limits on the neutrinoless double beta 
decay amplitude, we must have M rf - C ,M^ > 1 TeV for gQ e w ~ 0.2 [15]. This however does 
not constrain the slepton masses which could still be in the 100 GeV range. Unlike the 
conventional light neutrino mass and explicit R-parity violating contributions, these new 
contributions lead to RH polarization for the electron produced in the decay. 



Figure 13. New contribution to neutrinoless double beta decay due to — e c mixing. 
5.2 7T° — > e + e~ decay 

This new R-parity violating interaction also has interesting consequences for rare leptonic 
decays neutral pion and Kaon decays. We see from Eq. (2.5) that via t-channel u c exchange 
this leads to the process 7r° — > e + e~ with an amplitude given by A ~ g 2 9l w /M^ c . The 
current PDG bound [38] on this process is Br(7r° — > e + e~) < 6 x 10~ 8 . Using the bounds from 
neutrinoless double beta decay, we predict that in our model we have Br(vr° -)• e+e") < 10~ 8 . 
Note that if there is mixing in the right-handed charged current of the same order as the 
CKM mixings, then we would predict for the K — > e + e~ branching ratio at the level about 25 
times smaller than corresponding pion decay. This is about 3 times smaller than the current 
PDG quoted bound. In the LHC search described above, we already restrict ourselves to this 
allowed parameter range. 

6 Conclusion 

In summary, we have studied the phenomenology of a class of minimal SUSYLR models with 
dynamical R-parity breaking, i.e., R-parity must necessarily break in order for parity and gauge 




(5.2) 
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symmetry breaking to occur. This induces a new class of R-parity violating interactions due 
to the mixing between e c and w£, which are not present in the usual MSSM with explicit or 
spontaneous R-parity violation. These interactions lead to a new contribution to neutrinoless 
double beta decay which restricts the squark/gluino masses to be in the TeV range. The model 
has its characteristic signature at LHC which consists of final states of type e + e~ ^ [i jj or 
fj, e^bbjj for smuon as the NLSP. We estimate the background for this process and find that 
for M\y R not far above a TeV, the model should be testable once LHC reaches its full energy 
and luminosity. Incidentally, in this model there is also an upper limit on the mass of the 
right-handed W R boson in the low TeV range for symmetry breaking to occur. A large part 
of the mass range could be accessible even in the early running at the LHC. 
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A Fully parity symmetric version 

In this appendix, we consider the full parity symmetric version of the model. We now keep 
the A and A multiplets in our model of Table 1 . The Yukawa superpotential is given for this 
case by the same expression as Eq. (3.1) with two additional terms: L T T2AL and /iATrAA. 
The full potential that is parity symmetric is given below: 



The D-term potential as well as the scalar potential can be found in Refs. [2, 9]. The arguments 
for the existence of the dynamical R-parity breaking is same as in the parity asymmetric 




+ m\ Tr(A f A) + Tr(A ct A c ) + m\ Tr(A f A) + Tr(A ct A c ) 
+ i {M 2L X a L X a L + M 2R X a R X a R + M x X bl Xbl + M^X g X g ) 
+ Q t t 2 A\^t 2 Q c + Z T r 2 ,4f^T 2 L c + iA f (ZIt 2 A l L l + lfr 2 A R l R 
+ £ $afe Tr (r 2 ^r 2 6 ) + B A Tr (AA + A C A C ) + h.c. . 



) 



(A.l) 



version discussed in sec. 2. So we do not repeat this discussion here. The only question we 
address here is the status of a possible parity symmetric vacuum 3 with dynamical R-parity 
breaking. 

First, we would like to understand why the symmetry breaking in the SUSYLR model 
without Higgs triplets [8] is not compatible with the parity symmetry. The point is the 
LH and RH sneutrinos have opposite B — L charges, so the D-term potential contributes a 
negative cross term 



V D 



1 



c\2 



(A.2) 



which tends to minimize the potential in the parity conserving (y) = (y c ). This is why the 
authors of Ref. [8] have to start with parity asymmetric soft mass squared for sneutrinos. 
In contrast, the corresponding term in model with Higgs triplets becomes 



Vd 



j9bl 



2v\ + 24 



(u) 2 - 2v\ + 2v\ 



(A.3) 



where (A ) = vl and (A ) = vl- According the D-flat condition found out in Fig. 2, each 
bracket is very close to vanishing. Therefore, such D-term potential does not play significant 
role in forcing the vacuum to preserve parity and it is still possible to start with a symmetric 
potential. It has been was shown in [2] that if leptonic Yukawa couplings satisfy the bound 



Y, > 



2f{M\ - B A ) 



Ml 



(A.4) 



the parity violating minimum is indeed lower than the parity conserving. By choosing Ma, 
and B/\ appropriately, we can satisfy this bound so that the parity violating and R-parity 
violating minimum is the global minimum. 

B Explicit form of charged fermion mass matrix 



In this appendix, we present the explicit form of charged fermion mass matrix in the SUSYLR 
model. The spontaneous R-parity violation induces a mixing between the new chargino Wr, 
higgsino A c and the usual electron field. 

To see this explicitly, first note that parity violation at the TeV scale requires spontaneous 
R-parity breaking at a similar scale, i.e., (z?g) ~ vr ~ vr. We can write down the charged 

fermion mass |\P T Mg\l/ + h.c, in the basis of ^ = [(Wr, A c+ , e c+ ), (W^, A c ~, e 



IT 



M 5 



M 
M T 



M 



M 



1/2 



-V2g R v R 



VzgRVR -jiA 



(B.l) 



We thank S. Spinner for raising this point. 
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Following the similar arguments below Eq. (2.2), one finds the physical electron field mass 
term can be written as 

C m = -em e (6 ee e c + 9 eW W+ + 9 eA A c+ ) + h.c. 

= -em e e c + h.c. , (B.2) 

where 9 ee ,9 e w,0 e A are order 1 mixing parameters. 

In this model, the role played by A c , A c Higgses is to give mass to the RH neutrinos. 
Meanwhile, their superpartners enter in the above mixing matrix, but it does not change the 
generic prediction of large e c —W£ mixing. 
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